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1.  Introduction 


Rapid  experimental  advances  in  the  field  of  quantum  information  processing  make 
it  essential  to  understand  the  methods  for  verifying  and  benchmarking  quantum 
states  and  systems.  The  process  of  characterizing  unknown  quantum  states,  termed 
quantum  state  tomography,  consists  of  a  series  of  repeated  measurements  in 
different  bases  on  an  ensemble  of  identical  states.1,2  The  statistical  properties  of 
these  measurement  results  can  then  be  used  to  reconstruct  the  form  of  the  unknown 
state  using  a  maximum  likelihood  technique.1,2  Experimental  realizations  of 
quantum  state  tomography  have  been  implemented  in  the  characterization  of 
2-photon  pure  and  mixed  states, 3-5  3-qubit  Greenberger-Horne-Zeilinger  (GHZ)6 
and  W  states,7  and  qudits  in  the  orbital  angular  momentum  setting.8,9 

Quantum  state  tomography  is  a  very  general  technique  that  can  be  used  to 
characterize  the  quantum  states  of  any  number  of  particles  and  dimensions. 
However,  this  technical  report  focuses  on  the  2-qubit  scenario,  since  it  corresponds 
to  the  quantum  states  generated  in  our  experimental  testbed.10 


The  working  principles  of  2-qubit  state  tomography  are  most  easily  understood 
when  density  matrices  are  expanded  into  the  Pauli  matrices 

3 


P  =  Sij 


(1) 


i,j 


where  the  sigma  are  the  Pauli  matrices  given  by 


<7n  = 


(2) 


The  Sij,  referred  to  as  the  multiple-qubit  Stokes  parameter,  can  be  expressed  in 
terms  of  a  combination  of  local  projective  measurements  on  each  qubit  as 
$i,j  —  where  each  single  qubit  Stokes  parameter  is  given  by 


S0  =  ((p3\p\(p3)  +  Ws  \p\<Pz), 


(3) 


Si  =  {(pi\p\(pi)  ~  {(pi  \p\<pi),  (4) 

^2  =  ((p2\p\(pz)  -  (<p2  lpl<P2>>  (5) 


and 
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^3  =  <<P3IPI<P3>  ~  (<P3  IPI<P3>- 


(6) 


The  quantum  states  being  projected  onto  are  given  by  the  eigenvectors  of  the 
corresponding  Pauli  matrix 


\<pi)  = +  \v», 

W2)  =  ^m+i\v)), 


W3)  =  \H),  (9) 

and  their  orthogonal  compliments.  The  expected  measurement  probability  of  each 
of  these  combinations  for  a  given  p  is 

Si,j  =  Tr[(di®dj)p\.  (10) 

It  is  clear  from  the  above  formalism  that  an  experimentally  determined  set  of  joint 
measurement  probabilities  can  be  used  to  define  an  unknown  quantum  state.  (For 
simplicity,  we  have  restricted  this  discussion  to  projectors  based  on  the  Pauli 
matrices,  but  more  general  state  sets  can  be  projected  onto,  including  those  which 
are  not  orthogonal.)1,2 

However,  the  experimental  determination  of  the  Sy  requires  a  large  number  of 
measurements  to  approximate  the  true  probabilities  and  thus  reconstruct  the  actual 
density  matrix.  Consequently,  this  naive  method  is  prone  to  errors  for  small  sample 
sizes,  and  can  even  result  in  the  construction  of  unphysical  density  matrices.  For 
this  reason,  a  more  sophisticated  approach  based  on  maximum  likelihood 
estimation  is  used  to  define  physical  density  matrices  that  are  most  likely  to  produce 
the  recorded  measurement  results.1,2 

The  purpose  of  this  work  is  to  determine  how  statistical  counting  errors  (which 
prohibit  the  exact  determination  of  measurement  probabilities)  propagate  through 
our  tomography  system,  including  the  maximum  likelihood  estimation.  Unlike 
previous  research  that  studies  general  expressions  for  error  propagation  in 
entanglement  metrics,2  we  focus  on  specific  experimental  scenarios  relevant  to  our 
quantum  networking  testbed.  To  this  end,  we  perform  quantum  state  tomography 
measurements  for  5  distinct  experimental  scenarios  and  digitally  add  uncorrelated 
noise  to  these  measurement  results.  We  then  compare  the  reconstructed  density 
matrices  found  before  and  after  the  addition  of  noise,  and  determine  how  this  noise 
translates  into  errors  in  common  entanglement  measures.  These  results  allow  us  to 
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establish  which  quantity  variations  calculated  from  reconstructed  density  matrices 
are  a  natural  result  of  statistical  noise  and  which  signify  other  sources  of 
experimental  error. 

2.  Methods 


We  performed  our  analysis  using  quantum  state  tomography  software  made  freely 
available  by  the  Kwiat  Quantum  Information  Research  Group.11  This  software 
consists  of  MATLAB  functions  capable  of  performing  maximum  likelihood 
tomography,  which  finds  the  best  estimate  of  a  quantum  state  given  a  particular 
data  set.  Specifically,  we  employed  the  simple_2q2d_tomography  function,  which 
performs  this  procedure  using  the  measurement  bases,  the  singles  counts  at  each 
basis,  the  coincidence  counts,  and  the  accidental  counts.  We  also  used  several 
standard  entanglement  measures  implemented  in  this  software  through  the  function 
general_bell^fidelity,  which  finds  the  fidelity  of  an  input  density  matrix  with  its 
closest  pure  maximally  entangled  state,  as  well  as  concurrence ,  tangle,  and 
linear jentropy,  which  calculate  the  measures  indicated  by  their  names. 

The  singles  counts  (SI  and  S2)  are  the  number  of  events  for  which  one  detector  (D1 
or  D2)  “clicks”  (i.e.,  detects  a  photon).  Given  a  photon  pair  source  (such  as  our 
entangled  photon  source  (EPS)  shown  in  Fig.  1),  entangled  photons  and  noise 
photons  travel  in  both  channels  and  are  registered  as  singles  counts  when  detected. 
Coincidence  counts  (Cl 2)  occur  when  D1  and  D2  detect  photons  at  the  same  time. 
Since  the  EPS  creates  photon  pairs  in  both  channels  simultaneously,  coincidence 
counts  occur  when  photons  generated  by  the  EPS  are  detected  by  both  D1  and  D2. 
A  portion  of  the  coincidence  counts  result  when  both  detectors  happen  to  detect 
uncorrelated  photons  at  the  same  time.  These  are  known  as  accidental  counts  (A  12) 
and  measure  the  occurrences  where  two  independent  photon  streams  happen  to 
arrive  at  detectors  D1  and  D2  simultaneously. 

To  understand  how  errors  propagate  through  the  relevant  functions,  we 
systematically  added  noise  to  different  data  sets  and  calculated  its  affect  on  the 
output  density  matrices  and  accompanying  entanglement  metrics.  Noise  was 
inserted  by  adding  a  random  number  of  counts  to  the  singles  (SI,  S2)  and 
coincidence  (Cl 2)  counts  of  the  input  data.  These  random  counts  were  sampled 
from  a  normal  distribution  with  a  mean  of  zero  and  a  standard  deviation  equal  to 
the  square  root  of  each  individual  count.  Accidental  counts  were  then  calculated 
based  on  the  simple  relation  (A12  =  SI  x  S2  /  gates),  where  the  singles  counts 
include  the  random  noise  (gates  is  a  user-entered  setting  which  defines  how  long 
the  detectors  measure  the  EPS  output).10  These  values  were  then  used  in  the 


Approved  for  public  release;  distribution  is  unlimited. 

3 


tomography  functions  outlined  above.  Statistics  were  compiled  by  running  each 
numerical  test  100  times. 

Many  quantum  state  tomographies  were  performed  for  various  experimental  setups. 
The  primary  experimental  setups  that  were  tested  are  denoted  as  (N1)-(N4), 
including  (N2b),  for  a  total  of  5  distinct  cases.  Each  case  represents  a  different 
combination  of  pump  laser  power,  number  of  detection  gates,  applied 
polarization-dependent  loss  (PDL),  and  applied  polarization-mode  dispersion 
(PMD).  Figure  1  shows  schematic  diagrams  for  each  experimental  setup. 
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PDLE:  PDL  emulator 

CC:  Coincidence  electronics 


Fig.  1  Experimental  setup  for  each  of  the  five  tested  cases.  Case  (Nl)  depicts  a 
back-to-back  setup  with  low  pump  power.  Cases  (N2)  and  (N2b)  depict  the  setup  to  perform 
compensation  of  PDL  in  one  channel  by  applying  PDL  in  the  other  channel.  Case  (N3)  also 
depicts  compensation  of  PDL;  however,  channel  A  includes  additional  PMD.  Case  (N4)  depicts 
a  back-to-back  setup  with  high  pump  power. 

The  (Nl)  case  represents  a  “back-to-back”  measurement  with  a  low  pump  power 
setting  of  3100  in  the  Entangled  Photon  Analyzer  (EPA)  software  (resulting  in 
p  ~  .001)  and  a  long  measurement  time  of  300  million  detection  gates.  A 
back-to-back  measurement  is  one  in  which  a  transmitter  (such  as  our  EPS)  is 
directly  connected  to  a  receiver/measurement  apparatus  (such  as  our  polar 
analyzers  and  correlated  photon  detection  system)  to  characterize  the  behavior  of 
the  system  when  it  is  not  influenced  by  additional  transmission  effects.  As  such, 
the  back-to-back  measurement  denoted  by  (Nl)  has  neither  applied  PDL  nor  PMD. 
The  (N2)  case  represents  the  application  of  approximately  5  dB  PDL  to  both 
channels  when  using  a  high  pump  power  setting  of  3500  in  the  EPA  software 
(resulting  in  p  ~  .05)  and  a  short  measurement  time  of  30  million  detection  gates. 
The  (N2b)  case  differs  from  (N2)  by  only  one  parameter — a  detection  time  of  100 
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million  gates.  (These  2  cases  represent  experimental  demonstrations  of  PDL 
compensation.)  The  (N3)  case  differs  from  (N2b)  in  that  the  PDL  emulator  in 
channel  1  is  replaced  with  a  different  PDL  emulator  that  also  introduces  6.6  ps  of 
PMD.  (This  case  represents  experimental  demonstration  of  PDL  compensation 
when  subjected  to  decoherence  due  to  PMD.)  Finally,  (N4)  corresponds  to  a 
back-to-back  measurement  with  high  power  (as  with  (N2),  (N2b),  and  (N3))  and 
long  measurement  time  (300  million  gates,  as  with  (Nl)). 

The  goal  of  the  tomography  analysis  is  to  determine  how  statistical  fluctuations  in 
the  measured  number  of  counts  affect  the  calculated  density  matrix  and  its 
concurrence  and  purity.  Specifically,  calculation  of  the  density  matrix  relies  on  the 
difference  between  the  number  of  coincidence  counts  (for  two  photons  aligned  in 
the  same  polarization  basis)  and  the  accidental  counts.  This  value  was  therefore 
analyzed  for  each  of  the  experimental  setups  (C12eff  =  C12max  +  C12min  -  2  x  A 12). 
Table  1  includes  all  of  the  relevant  experimental  parameters,  measured  counts,  and 
calculated  results. 

Table  1  Relevant  experimental  parameters  (green),  measured  quantities  (orange),  and 
calculated  results  (blue).  The  calculated  results  include  the  effective  coincidences  and  the  noise 
of  relevant  density  matrix  elements  and  entanglement  metrics. 


h 

#  gates 

Si 

S2 

C12max 

C12min 

A12 

(Nl) 

0.001 

3.00E+08 

9.94E+04 

7.31E+04 

442 

39 

29 

(N2) 

0.050 

3.00E+07 

3.53E+04 

4.72E+04 

225 

30 

56 

(N2b) 

0.050 

1.00E+08 

1.15E+05 

1.55E+05 

786 

86 

184 

(N3) 

0.050 

1.00E+08 

8.31E+04 

1.08E+05 

407 

53 

96 

(N4) 

0.050 

3.00E+08 

1.59E+06 

9.47E+05 

31284 

5585 

5288 

C12efr 

(T(pil) 

tr(p44) 

a(pi4) 

"(C) 

"(P) 

(Nl) 

423 

0.0082 

0.0089 

0.0109 

0.0230 

0.0223 

(N2) 

143 

0.0177 

0.0108 

0.0173 

0.0517 

0.0416 

(N2b) 

504 

0.0088 

0.0102 

0.0149 

0.0202 

0.0179 

(N3) 

268 

0.0140 

0.0140 

0.0212 

0.0389 

0.0306 

(N4) 

26293 

0.0013 

0.0012 

0.0016 

0.0032 

0.0034 

3.  Verification  of  Data  and  Results 


To  verify  that  all  measured  coincidence  and  accidental  counts  agree  with  their 
expected  values,  a  quick  loss  budget  calculation  was  performed  to  compare  the 
relative  values  of  SI,  S2,  C12,  and  A12.  Since  A12  is  computed  directly  from  the 
singles  counts  (SI  and  S2),  the  measured  value  of  A12  was  verified  with  respect  to 
the  corresponding  S 1  and  S2  values.  The  calculated  values  of  these  parameters  were 
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found  to  be  in  good  agreement  with  the  measured  values  for  all  5  experimental 
cases,  as  discussed  below. 


First,  we  compared  the  values  for  all  of  the  high-power  cases  ((N2)-(N4)). 
Beginning  with  the  singles  counts  for  (N4),  analysis  of  the  loss  budget  showed 
agreement  between  the  singles  and  accidental  counts  for  cases  (N2),  (N2b),  and 
(N3).  To  compute  the  expected  values  for  (N3),  we  considered  the  change  in  counts 
due  to  the  reduction  of  the  detection  time  from  300  million  gates  (for  (N4))  to  100 
million  gates.  This  reduction  in  measurement  time  by  a  factor  of  3  decreases  the 
number  of  singles  counts  by  a  factor  of  3  (equivalent  to  -4.8  dB  of  loss).  Next,  we 
accounted  for  the  addition  of  5  dB  of  PDL  in  channels  1  and  2,  resulting  in  an 
additional  loss  of  -2.5  dB  in  both  S 1  and  S2.  Finally,  we  accounted  for  the  insertion 
loss  of  the  additional  equipment  in  each  channel  (the  polarization  controllers)  and 
the  addition  of  PMD  in  channel  1.  The  combined  loss  due  to  PMD  and  the 
polarization  controller  in  channel  1  was  measured  and  found  to  be  -4.8  dB. 
Therefore,  SI  for  the  (N3)  case  should  differ  from  the  (N4)  value  of  1.59E6  by 
-11.8  dB.  This  results  in  an  SI  of  1.05E5,  which  is  in  good  agreement  with  the 
measured  value  of  8.3 1E4.  Channel  2  does  not  include  any  added  PMD,  and  the 
insertion  loss  of  the  polarization  controller  was  measured  to  be  -1 .25  dB .  Therefore, 
S2  should  differ  from  the  (N4)  value  of  9.47E5  by  -8.55  dB.  This  results  in  a  final 
expected  value  of  S2  for  the  (N3)  case  of  1.33E5,  which  is  also  in  agreement  with 
the  measured  value  of  1.08E5.  Since  accidental  counts  are  due  to  simultaneous 
occurences  of  2  streams  of  singles  counts  (A12  =  SI  x  S2  /  gates),  the  agreement 
of  the  expected  singles  counts  with  the  measured  counts  reinforces  agreement  of 
the  accidental  counts. 

Next,  the  expected  values  for  the  (N2)  case  were  calculated.  We  first  considered 
the  change  in  counts  due  to  the  reduction  of  the  detection  time  from 
300  million  to  30  million  gates.  This  reduction  in  detection  time  decreases  the 
number  of  singles  counts  by  a  factor  of  10  (equivalent  to  -10  dB  of  loss).  We  then 
account  for  the  addition  of  5  dB  of  PDL  in  channels  1  and  2.  This  results  in 
increased  loss  of  -2.5  dB  in  channels  1  and  2.  Finally,  we  account  for  the  insertion 
loss  of  the  polarization  controllers  (there  is  no  PMD  for  this  case).  The  loss  of  the 
polarization  controller  in  channel  1  is  -3  dB,  resulting  in  a  total  decrease  in  SI  of 
-15.5  dB.  Therefore,  the  final  expected  value  of  SI  for  the  (N2)  case  is  4.45E4, 
which  is  in  agreement  with  the  measured  value  of  3.53E4.  The  insertion  loss  of  the 
polarization  controller  in  channel  2  is  -1.25  dB,  resulting  in  a  total  decrease  of 
-13.75  dB  in  S2.  Therefore,  the  final  expected  value  of  S2  for  the  (N2)  case  is 
3.98E4,  which  is  also  in  reasonable  agreement  with  the  measured  value  of  4.72E4. 
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The  (N2b)  case  analysis  differs  slightly  from  the  (N2)  case — the  calculated  SI  and 
S2  values  for  the  (N2)  case  are  multiplied  by  3.33  (to  account  for  the  increase  from 
30  million  gates  ((N2))  to  100  million  gates  ((N2b)).  This  results  in  calculated 
values  of  SI  =  1.48E5  and  S2  =  1.33E5,  which  are  in  agreement  with  the  measured 
values  of  1.1 5E5  and  1.55E5,  respectively.  Comparison  of  the  (N2)  and  (N2b)  cases 
shows  that  all  SI,  S2,  C12,  and  A12  values  differ  only  by  a  factor  of  3.33  because 
of  the  change  in  gate  number. 

Because  the  number  of  generated  photons  (and  therefore  the  values  of  S 1,  S2,  C12, 
and  A 12)  is  a  complicated  function  of  power  due  to  the  presence  of  multiple 
photon-generating  effects12-15,  a  “back  of  the  envelope”  calculation  is  inadequate 
for  rapid  verification  of  the  relationship  between  the  low-power  (Nl)  and 
high-power  ((N2)-(N4))  measured  counts.  However,  Jones  et  al. 12,13  show  that  the 
relative  number  of  counts  between  the  low  and  high  power  cases  is  in  agreement 
with  the  expected  values.  Similarly,  this  presence  of  multiple 
photon-generating  effects,  as  well  as  the  difference  in  applied  PDL  and  PMD16-20 
among  the  5  cases,  makes  loss  budget  calculations  such  as  those  performed  above 
insufficient  to  quantitatively  analyze  the  relative  magnitude  of  C12  for  each  of  the 
5  cases.  Since  the  coincidence-to-accident  ratio  (CAR)  (Cl 2/A  12)  provides  a  good 
measure  of  the  correlation  between  two  photon  streams,  we  use  it  to  verify  the 
relative  magnitude  of  C12.  For  example,  2  completely  uncorrelated  photon  streams 
will  have  a  CAR  of  1  because  all  measured  coincidence  counts  (Cl 2)  are  merely  a 
result  of  the  statistical  coincidence  that  1  photon  in  each  channel  is  detected  at  the 
same  time  (i.e.,  all  measured  coincidences  are  due  to  accidental  coincidences 
(A12)).  A  CAR  greater  than  1  means  that  the  photons  are  correlated  such  that  they 
are  being  generated  simultaneously,  and  more  than  the  expected  number  of 
coincidences  are  being  measured  for  2  independent,  uncorrelated  photon  streams. 

The  (Nl)  case  has  the  greatest  CAR  (13.9),  followed  by  the  (N4)  case  (5.3),  then 
the  (N2),  (N2b),  and  (N3)  cases  (~2.5)).  As  expected,  the  (Nl)  case  has  the  highest 
CAR  because  it  is  the  only  case  with  a  low  pump  power.  As  the  pump  power 
increases,  more  uncorrelated  photons  are  generated  in  channels  1  and  2  due  to 
Raman  scattering12,13'21-24  and  the  relative  amount  of  measured  coincidences  due  to 
entangled  photon  pairs  decreases.  The  (N4)  case  has  the  second  highest  CAR  due 
to  the  lack  of  applied  PDL  or  PMD.  Additional  PDL  and/or  PMD  is  applied  to 
channels  1  and  2  in  the  (N2),  (N2b),  and  (N3)  cases,  resulting  in  decoherence  of  the 
entangled  state  generated  by  the  EPS.  This  decreases  the  amount  of  coincidences 
measured  due  to  entangled  pairs,  which  reduces  the  CAR.  The  loss  budget 
calculations  used  to  determine  SI,  S2,  and  A12  and  the  qualitative  description  of 
the  CAR  verify  the  relative  magnitudes  of  all  of  the  SI,  S2,  C12,  and  A12  values 
for  each  of  the  5  cases. 
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Once  all  of  the  measured  values  (shown  in  the  orange  columns  in  Table  1)  were 
verified,  we  input  the  experimental  data  to  the  quantum  state  tomography  software 
discussed  in  Section  2  (the  results  are  shown  in  the  blue  columns  of  Table  1).  The 
relative  results  of  each  case  agree  with  the  general  trends  that  would  be  expected 
with  the  addition  of  random  statistical  noise;  that  is,  the  cases  with  a  greater  number 
of  measured  counts  have  lower  noise  in  the  calculated  density  matrix  elements 
(pi i,  p44,  p  14) ,  concurrence  (C),  and  purity  (P).  The  very  small  amount  of  statistical 
noise  in  all  of  these  parameters  confirms  that  statistical  noise  is  not  the  dominant 
cause  of  variation  in  equipment  performance  in  our  quantum  networking  testbed. 
For  example,  we  find  no  greater  than  approximately  0.05  standard  deviation  in  the 
concurrence,  but  we  often  perform  consecutive  quantum  state  tomographies  and 
find  a  standard  deviation  of  greater  than  0.05  in  the  concurrence.  We  postulate  that 
the  greater  variation  in  our  experimental  results  is  largely  due  to  drifts  in  the 
behavior  of  the  source,  including  power  drifts  and  temperature  drifts  (which  change 
the  polarization  characteristics  of  our  all-fiber  system). 

4.  Conclusion 


Quantum  state  tomography  is  a  powerful  tool  for  characterizing  and  benchmarking 
quantum  networking  systems.  This  technical  report  describes  the  effects  of 
statistical  counting  noise  on  density  matrix  elements  and  entanglement  measures 
derived  from  states  reconstructed  from  maximum-likelihood  tomography. 
Specifically,  we  have  focused  on  states  corresponding  to  different  experimentally 
relevent  scenarios  found  in  our  quantum  networking  testbed.  Artificially  adding 
uncorrelated  noise  proportional  to  measured  counts  did  not  cause  the  density  matrix 
elements  to  vary  with  a  standard  deviation  of  greater  than  approximately  0.05. 
Additionally,  no  standard  deviations  in  the  concurrence  and  purity  (derived  from 
the  density  matrices)  greater  than  approximately  0.05  were  found  in  each  case  and 
in  some  cases,  were  as  low  as  approximately  0.003.  From  these  results,  we  conclude 
that  statistical  counting  noise  is  not  a  major  contributor  to  variations  between  runs 
in  our  quantum  networking  testbed.  This  information  will  help  us  analyze  system 
performance  as  we  begin  to  implement  different  quantum  information  processing 
tasks25,26  using  our  quantum  networking  testbed. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


ARL 

US  Army  Research  Laboratory 

CAR 

coincidence-to-accidental  ratio 

EPA 

Entangled  Photon  Analyzer 

EPS 

entangled  photon  source 

PDL 

polarization-dependent  loss 

PMD 

polarization-mode  dispersion 
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